Transport Coefficients of Hadronic Matter near Tc 
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A hadron resonance gas model including all known particles and resonances with masses m < 2 
GeV and an exponentially rising density of Hagedorn states for m > 2 GeV is used to obtain an 
upper bound on the shear viscosity to entropy density ratio, rj/s ~ l/(47r), of hadronic matter near 
Tc- We found a large trace anomaly and small speed of sound near Tc, which agree well with recent 
lattice calculations. We comment on the bulk viscosity to entropy density ratio close to Tc- 
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The large azimuthal asymmetry of low-pT particles 
and the strong quenching of high-pT probes measured 
at RHIC [H indicate that the new state of matter pro- 
duced in heavy ion collisions is a strongly interacting 
quark-gluon plasma The matter formed in these 

collisions behaves almost as a perfect liquid [l] charac- 
terized by a very small value for its shear viscosity to 
entropy density ratio, which is in the ballpark of the 
lower bound rj/s > l/(47r) [4] derived within the anti- 
de Sitter/conformal field theory (AdS/CFT) correspon- 
dence Q. It was further conjectured by Kovtun, Son, 
and Starinets (KSS) that this bound holds for all sub- 
stances in nature (see, however, Refs. 0, [1| for possible 
counterexamples involving nonrelativistic systems). 

Recent lattice calculations ^ in pure glue SU (3) gauge 
theory have shown that rj/s remains close to the KSS 
bound at temperatures not much larger than T^. Addi- 
tionally, calculations within the BAMPS parton cascade 
[lo| . which includes inelastic gluonic gg <-!■ ggg reactions, 
showed that ry/s ^ 0.13 in a pure gluon gas fllj. More- 
over, it was argued in [l3 | that this ratio should have a 
minimum at (or near) the phase transition in quantum 
chromodynamics (QCD). This is expected because rjjs 
increases with decreasing T in the hadronic phase [l3[ 
(because the relevant hadronic cross section decreases 
with T) while asymptotic freedom dictates that 77/s in- 
creases with T in the deconfined phase since in this case 
the coupling between the quarks and the giuons (and the 
transport cross section) descreases logarithmically [11]. 
Note, however, that in general perturbative calculations 

(see. 



15|) 



are not reliable close to Tc (see, however, Ref. 

Thus far, there have been several attempts to compute 
77/s in the hadronic phase using hadrons and resonances 
. However, these studies have not explicitly 



considered that the hadronic density of states in QCD 
is expected to be ~ cxp(rn/Tf/) for sufficiently large m 
[H where Th ^ 150 - 200 MeV is the Hagedorn 
temperature [3| (see |2l| for an update on the exper- 
imental verification of this asymptotic behavior). This 
hypothesis was originally devised to explain the fact that 
an increase in energy in pp and pp collisions does not 
lead to an increase in the average momentum per parti- 



cle but rather to production of more particles of different 
species [l9[. Moreover, hadron resonance models that in- 
clude such rapidly increasing density of states are known 
to have a "limiting" temperature, Tm. nx, beyond which 
ordinary hadronic matter cannot exist [l9| . 

In this letter, a hadron resonance gas model that in- 
cludes all known particles and resonances with masses 
TO < 2 GeV 24I and also an exponentially increasing 
number of Hagedorn states (HS) [2^, [i^l is used to pro- 
vide an upper limit on 77/s for hadronic matter close to Tc 
that is comparable to l/47r. Additionally, we show that 
our model provides a good description of the recent lat- 
tice results [25] for the trace anomaly and also the speed 
of sound, Cs, close to Tc = 196 MeV. We comment on the 
effects of including HS on the bulk viscosity to entropy 
density ratio, C/s, of hadronic matter near Tc- 

The assumption behind hadron resonance models is the 
description of thermodynamic properties of a hadronic 
interacting gas by a free gas with these hadrons and 
their respective resonances. In [2^ the pressure of an 
interacting pion gas calculated within the virial expan- 
sion nearly coincides with that of a free gas of pions 
and p mesons. There is nearly exact cancellation be- 
tween the attractive and repulsive S-wave channels 27]. 
We assume that attractive interactions can be described 
by the inclusion of resonances which for large masses 
follow a Hagedor n sp ectrum. The system's mass spec- 
trum is assumed to be a sum over discrete and 
continuous states p{m) = phg("^) + PHsinT-), where 
Phg{^) = Sf^" gi^{^ ~ TOi) 6*(Mo — to) involves a sum 
over all the known hadrons |22| and their respective de- 
generacy up to Mq < 2 GeV [31] and for larger masses 
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where we take toq = 0.5 GeV, A = 0.5 GeV^ [231, and 
Th = Tc- In general, repulsive interactions among the 
hadrons soften the dependence of the pressure on the 
temperature 271, l28|, l29| . Their effects are included us- 
ing the excluded-volume approach derived in [j? ' 



where 

the volume excluded by a hadron equals its energy di- 
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FIG. 1: Comparison between e{T)/T* = (e - 3P) /T^ using 
our hadron resonance gas model including Hagedorn states 
with 2 < m < 20 GeV [S^l (solid red line) and our hadron 
gas model with only known hadrons up to m < 2 GeV (black 
dashed line). The blue band between the curves illustrates 
the effects of HS. Repulsive interactions are included via an 
excluded volume approach [2!i| with B^^'^ = 0.34 GeV. Lattice 
data for the p4 action with A^'^ = 6 25] is also shown. 



vided by AB, where B plays the role of an effective MIT 
bag constant. The thermodynamic quantities were found 
using 
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and the standard thermodynamic identities at zero 
baryon chemical potential [29|. Note that the temper- 
ature T and the pressure P{T) of the system (after vol- 
ume corrections) are defined in terms of the quantities 
computed in the point particle (subscript pt) approxima- 
tion (i.e., no volume corrections). When Ppt(Tc)/4:B < 1 
there is still a limiting temperature that is larger than 
Tc^. We take B^/* = 0.34 GeV in our calculations, 
which implies that T^ax > Tc- We restrict our discussion 
to T < Tc because at higher temperatures a description 
involving quarks and gluons should be more adequate. 

Our results for the trace anomaly are shown in Fig. [T] 
where the mass of the heavier Hagedorn state was set to 
be Mjnax = 20 GeV. Note that the inclusion of HS cor- 
rectly captures the trend displayed by the lattice data in 
the transition region whereas our hadron gas curve does 
not (s^l- This remains true if other values of B are used. 
We checked that our results did not change appreciably 
in this temperature range when Mmax is increased to 
80 GeV. This happens because the divergences normally 
associated with the limiting temperature only occur in 
this case at T^ax ~ 210 MeV. Were Th < Tc, the de- 
pendence of the thermodynamic quantities with Mmax 
would be much more pronounced. In general, a very 
rapid increase in the number of particle species (specifi- 
cally heavier species) around Tc is expected to strongly 
reduce the speed of sound = dP/ de at the phase tran- 
sition. While Cg — > at the transition would certainly 
lead to very interesting consequences for the evolution 
of the RHIC plasma |32], recent lattice simulations have 
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FIG. 2: ci including HS with 2 < m < 20 GeV (sohd red 
line), HS with 2 < m < 80 GeV (dotted-dashed blue), for a 
hadron gas model without HS (dashed black curve), and for 
p4 action lattice results with Nr = 6 [2^ (dotted curve). 



found that 



0.09 near It is shown in Fig. H 



that cl{T ~ Tc) - 0.09 in the model with HS while for 
the model without them ~ 0.25 near the transition. 
Note that when M^ax — 80 GeV (dashed blue curve) 

is only a bit smaller than 0.09 near the phase transi- 
tion. Other quantities such as the total entropy density 
near Tc are found to agree with lattice results within the 
uncertainties present in those calculations [sS]. 

The total shear viscosity of our multi-component sys- 
tem computed within kinetic theory [s^l is jytot ~ 
Q!^jni(pi)Ai, where rii is the number density, (pi) is 
the average momentum, and is the mean free path 
for discrete states and HS {a ^ Moreover, 

Xi — n.j aij^ where cry is the scattering cross sec- 
tion. Due to their large mass, the particle density of HS 
is much smaller than that of discrete states. Thus, one 
can neglect the small contribution to the mean free path 
from terms involving the interaction between the stan- 
dard hadrons and the HS. In this case, rjtot = Vhg + '>1hs 
where tihg is the shear viscosity computed using only 
the interactions between the standard hadrons while 
rjHs = 5 X)j ^i{p)i includes only the contribution from 
HS, which move non-relativistically since mns/T ^ 1. 
Note that our approximation for rjtot provides an upper 
bound since the inclusion of the interactions between HS 
and hadrons would decrease rjtot- Using the results above, 
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While the entropy dependent prefactor in Eq. ([3]) can be 
easily determined using our model, the detailed calcula- 
tion of r]HG and rjns requires the knowledge about the 
mean free paths of the different particles and resonances 
in the thermal medium. In the non-relativistic approx- 
imation, we can set (pi) ~ mi{vi) = y^Srrii T/n in Eq. 
([3]). Note that HS with very large m^'s are more likely 
to quickly decay. We assume that Xi — ti (vi) where 
T, = 1/r, = 1/(0.151 m, - 0.0583) GeV^i is the inverse 
of the decay width of the i*^ HS obtained from a linear 
fit to the decay widths of the known resonances in the 
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FIG. 3: 77/s is shown for a gas of tt's and nucleons fl^ (upper 
dashed black line), for a hadron resonance gas with excluded 
volume corrections [13 (lower dashed black line), and for KSS 
limit, r}/ s — l/4n, [g] (solid black). An upper bound from HS 
on r]/s (solid red) and the effects of HS (blue band) are shown. 



particle data book [23|, |2J, |35| . Our choice for gives 
the largest mean free path associated with a given state 
because it neglects any possible collisions that could oc- 
cur before it decays on its own. Note, however, that 
the decay cross section is in general different than the 
relevant collision cross section for momentum transport 
that contributes to rj according to kinetic theory. Thus, 
it is not guaranteed a priori that these decay processes 
contribute to i] in the usual way. Further studies of the 
relationship between HS and 77 could be done, perhaps, 
using the cross sections discussed in [36| . 

Substituting the results above in we find that rjns — 
8T^j niTi/Sn. The remaining ratio {ri/s)j^Q has been 



computed in Refs. [TtI . [TsI using different models and 
approximations. Since our main goal is to understand 
the effects of HS on {r]/s)tot, here we will simply use the 
values for {f]/ s)^q obtained in some of these calculations 
to illustrate the importance of HS. We chose to obtain 
{^1 s)hg for a gas of pions and nucleons from Fig. 5 in [l^ 
and for a hadron resonan ce g as with (constant) excluded 



volume corrections from [17|. Note that the results for 



•q/s obtained from the calculation that included many 
particles and resonances 17[ are already much smaller 
than those found in [l6j where only pions and nucleons 
are considered. A linear extrapolation of the results in 
1^ . 17 1 was used to obtain their 77/s values at high tem- 
peratures. In Fig. [21 (?7/s)j^( drops significantly around 
Tc because of HS. This result is especially interesting be- 
cause 77/s in the hadronic phase is generally thought to 
be a few times larger than the string theory bound. One 
can see that the contributions from HS should lower -q/s 
near to the KSS bound close to T^- Thus, the drop in 
77/ s due to HS could explain the low shear viscosity near 
Tc already in the hadronic phase. We used Mmax = 20 
GeV in Fig. [3] but the results do not change significantly 
if Mmax is increased by a factor of 4. 

The large value of the trace anomaly near Tc observed 
on the lattice has been used as an indication that C,/s oi 
QCD may be large at the phase transition 37, 3^. This 
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FIG. 4: Estimates for C/s = G^(0)/(9u;o s) (i^o = 1 GeV) for 
the model that includes HS with 2 < m < 20 GeV (solid red 
line) and 2 < m < 80 GeV (dotted-dashed blue line) and our 
hadron gas model with m < 2 GeV (black dashed line). 



is very different than at high temperatures where C/s is 
predicted to be small (39| . This may have some interest- 
ing phenomenological consequences such as the forma- 
tion of clusters at freeze-out 40]. Using QCD sum rules 
in [41J , one can extract the (zero- momentum) Euclidean 
correlator of the energy- momentum tensor's trace, 

G^(0,O) = J d^x {0^^iT,^)9;iO,O)) 

= (T5t - 4) (e - 3p) (4) 

According to [3^, C can be obtained via using 
TTp{uj, 0)/9uj = Cujq/{uj'^ + ujq) as an ansatz for the small 
frequency limit of the {90) spectral density at zero spa- 
tial momentum, p{u},0). The parameter ujo{T) defines 
the energy scale at which perturbation theory is applica- 
ble. The validity of this ansatz has been recently studied 
in Refs. [i^. Here we assume that this ansatz can at 
least capture the qualitative behavior of C around Tc and 
we use it to estimate how HS change the (/s close to Tc- 
The results for C/s = G^(0)/(9wos) are shown in Fig.g] 
where wq = 1 GeV. While C/s decreases near for the 
hadron gas model, C/s including HS increases close to Tc 
and this enhancement does not vary much with Mmax- 

In conclusion, a hadron resonance gas model includ- 
ing all the known particles and resonances with masses 
m < 2 GeV and also an exponentially rising level den- 
sity of Hagedorn states for m > 2 GeV was used to ob- 
tain an upper bound on rj/s for hadronic matter near Tc 
that is comparable to the KSS bound 1 / (47r) . The large 
trace anomaly and the small near Tc computed within 
this model agree well with recent lattice calculations [2^ . 
Moreover, according to the general result that small rj/s 
implies strong jet quenching [43| . our significant reduc- 
tion of rj/s indicates that hadronic matter near the phase 
transition is more opaque to jets than previously thought. 
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